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Abstract 


Detailed simulations of buffer gas cooling, which is widely applied at the 
rare-isotope facilities, have been investigated by using three frequently used 
models: viscous damping force (VDF) model, hard sphere collision (HSC) 
model, and realistic interaction potential (RIP) model. By comparing the 
stopping ranges from the simulated data with those from SRIM code, we 
found that the VDF and HSC models are better in low energy range, but the 
RIP model is better in relatively high energy range. In helium buffer gas, 
which is always used in cooler and bunchers, the VDF model should be used 
when the ion’s energy is less than ~5 eV/u, and the RIP model less than 
~80 eV/u. Above this energy, all three models can not be used. 
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1. Introduction 


Buffer gas cooling is widely applied at the rare-isotope facilities nowa- 
days [1-4]. By combining with radiofrequency ion guide, the continuous ion 
beams can be improved and the cooled ion bunches with low emittance and 
small spread can be provided. This method is universally applicable to all 
elements ranging from helium to the heaviest elements [5], and is rather fast 
that the ions could be cooled within single pass through the cooler device. 
It has become one of the essential parts for the successfully high-precision 
Penning trap mass measurements. 

In order to obtain the high performance of ion beam cooler, detailed sim- 
ulation of buffer gas cooling is critical. There are three models frequently 
used for such a simulation: viscous damping force (VDF) model, hard sphere 
collision (HSC) model, and realistic interaction potential (RIP) model. VD- 
F model is a macroscopic approach by providing a simple prescription of 
the time-averaged cooling force, and the others are microscopic methods by 
tracking the cooling process by simulating individual collision between the 
ion and the buffer gas molecule. 

The buffer gas cooling effect has been simulated by many researchers. 
Lunney and Moore [6] developed a simulation package to simulate a radio 
frequency quadrupole rod system operated in helium buffer gas for cooling ion 
beam to reduce its emittance. Results indicated that an emittance reduction 
of almost 2 orders of magnitude in all coordinates was possible with near 
100% transmission barring charge exchange and molecular formation. 

Kim [7] used both the HSC and RIP models to simulate the cooling 


process of ions in a gas-filled ion guide and to compare the calculated tem- 


perature of cooled ions to experimental data. Only the calculations with 
the realistic potentials were able to simulate both realistic cooling rates and 
reproduce the experimentally obtained temperature. 

Schwarz [8, 9] developed a Monte-Carlo simulation package to calculate 
the motion of ions in the presence of arbitrary electric (and/or magnetic) 
fields and buffer gas. The mobility data Ko as a function of drift velocity 
obtained from simulations with a realistic scatting potential were in good 
agreement with the experimental data [10, 11]. But this is not the case for 
the simulations with a hard-sphere collision. 

We emphasize that all these simulations have been performed for the 
ions’ energies less than 100-200 eV and not beyond. The suitable energy 
(velocity) range (or energy limitation) is still unknown for these models in 
the simulation of buffer gas. In this paper we would like to present our results 


by comparing the stopping ranges of various ions in different buffer gases. 


2. Models 


2.1. Hard-sphere collision (HSC) model 


The HSC model is the simplest approach to handle the cooling force, 
and it treats the collision between the ion and the buffer gas particle as a 
collision of hard spheres. The ion travels through the neutral background 
gas, and collides with those gas particles, which travel in randomized ve- 
locities according to the Maxwell-Boltzmann distribution, and then scatters 
elastically. This approach has been used in a number of calculations and 


implemented in several simulation packages [7, 12]. 


2.2. Viscous damping force (VDF) model 


Viscous damping is produced by the long-range induced-dipole interaction 
by which an ion feels the effect of many gas atoms at once rather than the 
very short range force of a direct collision. The experimental values of ion 
mobility at low energies enable more accurate simulation of viscous damping 
than hard-sphere collisions. The validity of VDF model has been verified by 
the cooling time of °°K in helium gas [6]. 

If an ion with a charge of q and a mass of m moves in gas, the ion motion 
is viscously damped, in principle down to the temperature of the gas itself. 
The time-averaged cooling effect of the buffer gas is often approximated by a 
viscous drag force. In the simulation, the effect of buffer gas cooling can be 
explained as ion’s acceleration adding a viscous damping acceleration term, 


which could be express as [6]: 


du? q du 1 Uth 
sd- 1 
di? m ( dt K T Kin , ( ) 


where u stands for the ion’s position, such as x, y and z. K is the ion’s 


mobility, and vn is the thermal velocity of the ion and Ky, its correspond- 
ing mobility. The last term is added to prevent the unrealistic situation in 
which the ions are cooled to zero temperature since the statistical nature of 
the collisions of ions with the buffer gas molecules has not been taken into 
account. 

The mobility K value at a given pressure of p and temperature of T can 


be normalized from the “reduced” value Ko by 


T 1013 mbar 


KER.. f 
0" 273.16K p 


(2) 


Table 1: Fitted parameters for VDF model for different ion/molecule combinations 
Ion/Mol a bo by b2 
Art/He 1.34x 10° 6.29x107 —8.21x10° 0.299 
Cst/He 2.41 x 104 2.22 x 108 —1.82x 10° 0.395 
Hg*/He 2.70 x 104 2.72x 10° -1.92 x 10° 0.360 


Cst/Ne 2.83x10° 1.09x 10” -5.12x 103° 1.043 
Cst/Ar 1.58x 10? 8.86x 10° -3.05x 10° 1.478 


Ko values are known for a wide variety of ions in the low energy part [10, 11, 


13]. For the high energy part, they can be approximately calculated from: 


2qg(m+M)? 1 
Ky = i 3 
2 m mM Nov’ (3) 


where M is the mass of the gas particle, N the gas density and o the cross 
section for hard-sphere collision. 


The whole curve can be easily fitted by the following formula: 


v+a 
Ko = f 4 
2 bo + bi (v + a) + bo(u + a)? (4) 


where a, bo, bı and bz are fitting parameters. As an example, the reduced 
mobility Ko of !°8C'st ion in helium gas is shown in Fig. 1. The data point 
with velocity lower than 10,000 m/s are taken from Ref. [13] and the others 
are calculated from Eq. (3). The fitted parameters for different ion-molecule 


combinations are listed in Table 1. 
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Figure 1: The reduced mobility Ko of 1°°Cs* in He-gas versus drift velocity. The data 
point with velocity lower than 10,000 m/s are taken from Ref. [13] and the others are 
calculated from Eq. (3). The fitting parameters used: a = 2.41 x 10%, bo = 2.22 x 108, 
bı = —1.82 x 104, bz = 0.395. 


2.8. Realistic interaction potentials (RIP) model 


Although the idea of collisions of hard sphere is quite useful to understand 
the cooling mechanism, it ignores the long-rang interaction between the ion 
and the buffer gas particle, which lets the gas act as a viscous medium and 
lets the ion experience a damping force. The RIP model! considers this 
interaction in a “realistic” way. 

The RIP model chooses the potentials V(r) (also called (n,6,4)-potentials) 
at a distance of r as the sum of one short-range repulsion term and two 


attractive terms [14], 


V(r) = —-—-—. (5) 


'The name “Realistic interaction potentials (RIP) model” may be not very appropriate, 


instead, we may call it n-6-4 model because the (n,6,4)-potentials are used. 


The C4/r* and Og/r® terms quantify the attractive interaction of the ion’s 
charge with the electric dipole and quadrupole moments, respectively. 


Eq. (5) can alteratively be written in the dimensionless form as 


ne 
n(3+7) —120+7) 


« [Porn (2) -a 


-3a -m ()'], 6) 


V(r) 


where rm and € are the position and depth of the minimum of the potential, 
respectively, and y is the dimensionless parameter that measures the strength 
of the r~® term relative to the r~* term. 

The parameters in the potentials (Eq. (5) or (6)) can be derived by fitting 
the so-called collision integral, a theoretical expression for the collision cross- 
sections, to experimental data obtained from mobility data [14, 15]. 


The mobility K can be represented by 


2 M 1 
K- 3q Tt m+ l (7) 
16N kTeff mM MTerr) 


where g and m are the charge and mass of the ion, respectively, M is the 


mass of the gas particle, N is the number density of the gas molecules, and 
k is Boltzmann’s constant. The effective temperature Teff accounts for both 


the gas temperature Tas and the drift velocity va of the ion: 


3 3 1 
ghlers = 5 Toas + 5M. (8) 


The collision integral (of first order) Q(T) is a function of temperature 


un = (Er 


x je — cos(d(b, E)))b dab dE. (9) 


The deflection angle 0(b, Æ) is calculated as a function of the impact param- 


eter b and the relative energy E of the collision: 


ai 1 d 
6(b, E) =a — 2» f - (10) 
Ta ze oe 
r2 E 
where the distance of closest approach ra is the outermost root of 
bd? Vra) 
Sas =0. 11 
es (11) 


By inserting Eq. (5) or (6) into Eq. (10) and the result into Eq. (7) one 
obtains a relationship between the sets of potential parameters (n, B,C6,C4) 
or (N, m,€,7) and the collision integral. 

The fitted parameters for the (n,6,4)-potentials for a number of ion- 
molecule combinations can be found in Ref. [9]. These potentials can also be 
compared with some high-level ab initio calculations. Larry A. Viehland et 
al. [16-20] have performed such calculations to investigate the interaction- 
s between the atomic cations and the rare gas particles. They calculated 
the potential energy curves and compared calculated spectroscopic param- 
eters and transport coefficients with the experimental mobility data. The 
agreements were excellent. Figure 2 shows the ab initio calculated potential 
energy curves for C'st/He and Hg*/Ar, along with the potentials calculated 
from Eq. (5) by using parameters in Ref. [9]. As may be seen, in general, 


the two sets of potential curves are in good agreement, and although the 
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Figure 2: The potential energy curves for Cst/He (upper panel) and Hg*/Ar (lower 
panel). The black dashed curves represent the potentials calculated from Eq. (5) by using 
parameters in Ref. [9], while the solid red curves represent the potentials taken from 


Ref. [21]. 


curves from the (n,6,4)-potentials change faster than that from the ab initio 
calculations at short ion-molecular separations. Among all the combinations 
we considered in this paper, the potential curves for Hg*/Ar combination 


have the best agreement. 


3. Simulations 


As the buffer gas cooling is typically used in the ion manipulation in the 
ion cooler and buncher, such as those in ISOLTRAP [22], SHIPTRAP [23], 
JYFLTRAP [24, 25], LEBIT [4], CPT [26], LPT [27, 28] and so on, not only 
the electric (and/or magnetic) fields should be known accurately, but also the 


effect of the presence of the buffer gas should be modeled adequately. In our 


simulations, after the electrode layout had been defined, the ion-optics sim- 
ulation package SIMION [29] was used to obtain the electromagnetic forces 
acting on the ions. However, for the purpose of this paper, we disabled, i.e. 
put 0 V on, all the electrodes, thus the movement of the ions is not affected 
by the presence of electrodes. 

To simulate the buffer gas cooling, we used three different models men- 
tioned above: VDF, HSC, and RIP models. In the scheme of the SIMION 
code, we developed a set of user programs that realized the VDF and HSC 
models to simulate collisions between the injected ion and the buffer gas. For 
the RIP model, it is very difficult to do the simulation by using the SIMION, 
thus we used the “IonCool” code [8], which was developed by Stefan Schwarz, 
instead. For convenience and comparison the “IonCool” was also used for 
the HSC model. After detailed simulations, we found that the HSC model 
results from both the SIMION and the “IonCool” were in good agreement, 
so we only reported here the results from the SIMION for the HSC model. 

For all models, the combinations of three species of ions, relatively light 
ion - * Ar*, medium ion - %C’st and heavy ion - ?°'Hgt, and three buffer 
gases, He, Ne and Ar, which are often used in the cooler and bunchers, had 
been simulated. The total number of ions for every simulation was 1,000. 
The traveling distance of the ion in buffer gas was chosen as the stopping 
range when the ion’s kinetic energy was decreased down to about its thermal 
energy. 

According to the study of Major and Dehmelt [30], only the results for 
those combinations that the ion’s mass is much heavier than the mass of 


buffer gas particle are reported here, since the cooling effect is our interest. 
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In our calculations, the parameters in Ref. [31], Ref. [9] and in Table 1 have 
been used for the HSC, RIP and VDF models, respectively. 

We used SRIM code [32] for the purpose of comparing the simulation re- 
sults with the experiments. By using their large collection of stopping power 
data, Helmut Paul and Andreas Schinner [33] compared those data to vari- 
ous stopping power tables and computer codes to estimate the reliability of 
these tables and codes. They found that SRIM worked fairly well for all ions 
and for all energies. Thus we chose the SRIM code to calculate the stopping 


powers/ranges and then used these calculated values as “experimental” data. 


4. Results and analysis 


The model differences are illustrated by Fig. 3, as an example, by com- 
paring calculated trajectories based on the VDF, HSC, and RIP models. The 
chamber is filled with helium buffer gas at a pressure of 1 Pa without the 
presence of any electric field, and one ?° Ar* ion is injected into the chamber 
along the Z axis with a kinetic energy of 400 eV. 

In the simulation by using the VDF model, the “Ar ion goes along the 
Z axis and there is no divergence to the X/Y axes (radial direction) except 
at the very end of the trajectory, and it almost stops at last. This is not 
the case for the HSC and RIP simulations. The ions go along the Z axis at 
first, then diverse to the radial direction, and at the end of the trajectories, 
the ions go randomly and never stop. The reason of this difference lies on 
whether the collision of ion with the gas particle has been taken into account, 
i.e., it is not taken into the VDF simulation, but it is in the HSC and RIP 


simulations. 
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Figure 3: Drift simulation of a single 4 Ar* ion in helium buffer gas at a pressure of 1 Pa 
by using different models. The ion is injected into the chamber, where no any electric field 
exists, along the Z axis with a kinetic energy of 400 eV. VDF denotes viscous damping 
force model, HSC hard-sphere collision model, and RIP realistic interaction potentials 


model. 


The corresponding kinetic energies of Ar ion versus elapsed time are 
shown in Fig. 4. It is obvious that the HSC simulation gets the shortest 
cooling time among all the three simulations, and the RIP the longest. The 
faster the energy losses, the shorter the cooling time and the smaller the 


stopping range. 


The stopping ranges of Art, Cs* and Hg* ions in He buffer gas obtained 
from simulations by using different models mentioned above are shown in the 
upper panels of Fig. 5. The results from SRIM are also shown for comparison. 


It is obvious that the VDF and HSC models are better than the RIP model 
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Figure 4: The kinetic energy of a “?Ar* ion in helium buffer gas versus elapsed time. See 


text and Fig. 3. 


when the ion’s energy is less than ~5 eV /u and the stopping ranges agree 
with the SRIM results better. In the low energy range, the RIP model 
underestimates the stopping power, i.e., overestimates the stopping range, 
of the ions in He buffer gas. But in the relatively higher energy range, the 
RIP model works better than the others. If we only compare the VDF model 
with the HSC model, we can find that the VDF model is much better. It can 
be explained by the experimental mobility values, which is only available for 
low energy ions, used in the VDF model. 

For the suitable energy range, generally speaking, in He buffer gas, the 
VDF model should be used when the ion’s energy is less than ~5 eV/u, and 
the RIP model less than ~80 eV/u. Above this energy, all three models can 


not be used and the errors from such simulations are too large. 


Figure 5 shows also the corresponding stopping ranges of Cs* ion in He-, 
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Figure 5: The stopping ranges of Ar* (upper left), Cs* (upper middle) and Hg* (upper 
right) ions in He-gas, and Cs* ions in Ne- (lower middle) and Ar-gas (lower right) as a 
function of ion’s energy per nucleon. The error of every data point is the statistical error 


from the simulation. For the VDF results the error bar is smaller than the symbol size. 


Ne-, and Ar-gas obtained from simulations. When the buffer gas becomes 
heavier, the upper limits to use the models decreases. For the RIP model, 
this value decreases from ~80 eV /u in He-gas, to ~50 eV/u in Ne-gas, and 
then to ~20 eV/u in Ar-gas. And for the VDF model, it decreases from ~5 
eV/u in He-gas, to ~2 eV/u in Ar-gas. 
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5. Conclusions 


The VDF, HSC and RIP models to simulate the behavior of ions in buffer 
gas have been investigated in detail. We found that 


e The VDF and HSC models are better in low energy range, but the RIP 


model is better in relative high energy range. 


e In He buffer gas, which is always used in cooler and bunchers nowadays, 
the VDF model should be used when the ion’s energy is less than ~5 
eV/u, and the RIP model less than ~80 eV/u. Above this energy, all 


three models can not be used. 


e If the buffer gas becomes heavier, the upper limits to use the models 
decrease. For the Cs* ion, the upper limit for the RIP model decreases 
from ~80 eV/u in He-gas, to ~50 eV/u in Ne-gas, and then to ~20 
eV/u in Ar-gas. And for the VDF model, it decreases from ~5 eV/u 


in He-gas, to ~2 eV/u in Ar-gas. 
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